Bulk viscosity of strange quark matter: Urea versus non-leptonic processes 
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A general formalism for calculating the bulk viscosity of strange quark matter is developed. 
Contrary to the common belief that the non-leptonic processes alone give the dominant contribution 
to the bulk viscosity, the inclusion of the Urea processes is shown to play an important role at 
intermediate densities when the characteristic r-mode oscillation frequencies are not too high. The 
interplay of non-leptonic and Urea processes is analyzed in detail. 
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I. INTRODUCTION 

Compact (neutron) stars provide a natural laboratory 
of matter under extreme conditions. In the central re- 
gions of such stars the baryon density of matter could 
reach values up to 10 times the nuclear saturation den- 
sity (i.e., 10po where po ~ 0.15 fm~ 3 ). At such high den- 
sity matter is likely to be in a deconfmed state in which 
quarks rather than hadrons are the natural dynamical 
degrees of freedom 0, 0, |, i, i] . 

One can argue that the ground state of deconfmed 
quark matter is a color superconductor. (For reviews 
on color superconductivity, see Refs. d @, II, H El EH, 
EH EH)- Many phases of color superconductivity are 
known that could possibly be realized in dense matter. 
It remains unclear, however, which of these describe the 
ground state of matter under the specific conditions in 
stars. This is because of theoretical uncertainties in treat- 
ing the strongly coupled, non-perturbative dynamics in 
QCD at the baryon densities of relevance. 

In order to clarify this one could use observational data 
from stars to narrow down the range of possibilities. For 
this program to work, one requires a detailed knowledge 
of the physical properties of various phases of matter that 
are likely to exist in stars. The transport properties are 
most suitable for developing sufficiently sensitive, as well 
as unique and verifiable predictions regarding observa- 
tional signals from stars. In this paper, we study the 
bulk viscosity of the normal phase of three-flavor quark 
matter. This work extends our recent analysis of the 
bulk viscosity in two-flavor quark matter [hij by includ- 
ing stran ge q uarks. This also augments existing studies of 
the topic HE El El El El M, El S3, in which the role 
of the Urea processes was not thoroughly investigated. 

In general, the bulk viscosity is a measure of the ki- 
netic energy dissipation during expansion and compres- 
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sion of a fluid. In compact stars, the density oscil- 
lations of interest have characteristic frequencies that 
are of the same order of magnitude as the stellar ro- 
tation frequencies. These are bound from below and 



from above, 1 s 1 < u> < 10 3 s 



(For the fastest- 



spinning pulsar currently known, PSR J1748-2446ad, one 



has i 



4.5 x 10 s corresponding to v = 716 Hz [23|.) 



The most important microscopic processes that provide 
the energy dissipation on the corresponding time scales 
are weak processes. Under conditions in stars, in partic- 
ular, the bulk viscosity of quark matter is determined by 
the combined effect of the flavor-changing weak processes 
diagrammatically shown in Fig. [TJ When an instanta- 
neous departure from chemical equilibrium is induced by 
expansion/compression of matter, the weak processes try 
to restore the equilibrium state and, while doing this, re- 
duce the oscillation energy. 

It is commonly argued that the bulk viscosity in the 
normal phase of three-flavor quark matter is dominated 
by the non-leptonic weak processes u + d <-» u + s 
[11 El ES El, 111 H EL EH . These are shown diagram- 
matically in Figs.[lja) and[TJb). As for the Urea (semi- 
leptonic) processes, see Figs.[]Tc)-QTf), they have consid- 
erably lower rates which are suppressed parametrically by 
a factor of order (T/ /x e ) 2 . This is in contrast to two-flavor 
(non-strange) quark matter, in which case the Urea pro- 
cesses u + e~ — > d + v and d — > u + e~ + v [see Figs. QTc) 
and QJf)] are the only ones that contribute [3]. (For 
studies of the viscosit y in various phases of dense nuclear 
matter, see Refs. 0, M HI EL EE H, M, M, H HI . ) 

In this paper, we study the bulk viscosity of three- 
flavor (strange) quark matter. One of the key issues that 
we address here is the interplay between the Urea and 
the non-leptonic processes. We shall show that, because 
a resonance-type phenomenon determines the bulk vis- 
cosity and because there is a subtle interference of the 
two weak processes, the simple argument about the dom- 
inance of the non-leptonic processes is not always justi- 
fied. 

This paper is organised as follows. In the next sec- 
tion, we develop the general formalism for calculating the 
bulk viscosity in three-flavor quark matter, paying spe- 
cial attention to the interplay of several different weak 
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5n = Suq Re(e lwt ) where Sno and u> are the magnitude 
and the frequency of the oscillations, respectively. The 
energy dissipation can be also expressed in terms of the 
induced pressure oscillations, 



(fdiss) = - / PVdt 
r Jo 



(2) 



where V = l/n is the specific volume, and r = 2tt/lu is 
the period of the oscillation. 

The density oscillations drive strange quark matter out 
of [3 equilibrium. While the weak processes tend to re- 
store equilibrium, they always lag behind. Thermal equi- 
librium, in contrast, is restored almost without any delay 
by very fast strong processes. The corresponding instan- 
taneous quasi-equilibrium state is unambiguously charac- 
terised by the total baryon density n, the lepton fraction 
X e , and the strangeness fraction X s , 



X, 

x s 



n 



+ n d + n s ) . 



(3a) 
(3b) 
(3c) 



where n u , rid, and n s are the number densities of up, 
down, and strange quarks, while n e is the number density 
of electrons. Because of the charge neutrality constraint, 
the number densities satisfy the following relation: 



1 



1 



-n d 



0. 



(4) 



FIG. 1: Diagrammatic representation of the j3 processes that 
contribute to the bulk viscosity of quark matter in stellar 
cores. 



processes. In Sec. [TTT] we calculate the bulk viscosity in 
the normal phase of strange quark matter. In the same 
section, we also study the role of the Urea processes on 
the bulk viscosity. The discussion of the results is given 
in Sec. HV1 



II. 



BULK VISCOSITY IN STRANGE QUARK 
MATTER 



In this section, we give a general derivation of the bulk 
viscosity resulting from the combined effect of all weak 
processes shown diagrammatically in Fig. [T] Following 
the same method as in Ref. [14| , one can relate the bulk 
viscosity £ to the average energy-density dissipation in 
pulsating matter, 



(£ diss )___l _ 



(1) 



(Note that there are no microscopic processes that would 
lead to deviations from this charge neutrality constraint.) 
Using this constraint together with the definitions in 
Eq. ([3]), one can express the separate densities and, in 
fact, any thermodynamic quantity of quasi-equilibrium 
quark matter in terms of n, X ei and X s . For the number 
densities, for example, one finds 



n e = X e n, 

n u = (1+X e )n, 

n d = (2 - X e - X s ) n. 

71, = A%77. 



(5a) 
(5b) 
(5c) 
(5d) 



The number densities can also be expressed in terms of 
the corresponding chemical potentials, n,; = ni(/ii). In )3 
equilibrium, the chemical potentials are related as fol- 
lows: ii s = iid = A*« + Me- I n pulsating matter, on 
the other hand, the instantaneous departure from equi- 
librium is described by the following two independent 
parameters: 



5^2 



Sfi s 

Me : 



- Spid, 



(6a) 
(6b) 



In the derivation of this relation, we assumed that small 
oscillations of the quark matter density are described by 



where S/ii denotes the deviation of the chemical poten- 
tial Hi from its value in (3 equilibrium. Note that the 
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third possible combination of the chemical potentials, i.e., 
8/13 = fid — Mit — He = S/u,2 — Shx, is not independent. 

In quasi-equilibrium, both quantities Shx and <5/x 2 can 
be expressed in terms of the variations of the three inde- 
pendent variables, 8n, 8X e , and 8X Sl 

8n 

Sfi, = d— + Bi5X e + Ai5X s , (i=l,2) (7) 
n 

where, as follows from the definition, the coefficient func- 
tions are given by 



A t 


\ dn d dn s J ' 




(8a) 


Bi 






(8b) 


C x 


d\x s dfi d 

= n s - lid- — , 

dn s an d 




(8c) 


A 2 


= n 7T~ ' 
an s 




(8d) 


B 2 


( dfX u dfl e \ 

\ dn u dn e ) 




(8e) 


c 2 


du s dfi u 
on s on u 


neT: — • 


(8f) 



When the quantity Sh% is non-zero, the two weak pro- 
cesses in the ith row of Fig.Q]have slightly different rates. 
To leading order in 8 Hi, we could write 



r (a ) 


-r (6 ) = 


-XxShx, 


(9a) 


r (c) 


-r (d) = 


-X 2 8h 2 , 


(9b) 


P (e) 


^ r (/) = 


-A 3 (6n 2 - 6 (a) , 


(9c) 



(Note that our conventions are such that the quantities 
Aj are non-negative.) The net effect of having different 
rates for each pair of the weak processes in Fig. [1] is a 
change of the electron and strangeness fractions in the 
system, i.e., 



d(8X e ) 

1 dt 
d(8X s ) 
1 dt 



A 2 <5/x 2 + A3 (Su2 — Su\) , (10a) 
-\ 1 Su 1 - A 2 <fyi 2 . (10b) 



According to these equations, the instantaneous val- 
ues of the electron and strangeness fractions, i.e., X e = 



Xe ' 1 + SX e and X s = Xg'' + 5X S , tend to approach their 
equilibrium values, xjP^ and xi°\ respectively. For ex- 
ample, a deficit of electrons (indicated by either Sfj. 2 > 
or S/j.2 — 8 fix > 0, or both) causes their production, see 
Eq. pOap . Similarly, a surplus of strange quarks (indi- 
cated by either Sfii > or Sfi 2 > 0, or both) is to be 
removed through the weak processes, see Eq. ()10bl) . Be- 
cause of finite rates for the weak processes, however, the 
oscillations of 8X e and 8X S always lag behind the density 
oscillations. In order to see this explicitly, we substitute 
8 \Xi from Eq. |7]) into Eq. (jTUJ) , and get the following set 
of equation for SX e and 8X S : 



d(8X e ) 

dt 
d(8X e ) 

1 dt 



Sn 

[(A a + A 3 ) C 2 - A3C1] — + [(A 2 + A 3 ) B 2 - X3B1] 8X e + [(A 2 + A 3 ) A 2 - X3A1] SX S , (11a) 
n 

8li 

- [Aid + A 2 C 2 ] [XiBt + \ 2 B 2 ] 8X e - [XiAx + X 2 A 2 ] SX S . (lib) 

n 



r 



The periodic solution to this equation can be found most 
easily by making use of complex variables. Thus, by 
defining 



SX e = Re (8X e , e lu}t ) , 



Re (8X S 



e 



(12a) 
(12b) 



we derive the following results for the complex magni- 
tudes: 



8X e0 = 



8X sQ = 



8tiq d\ + id 2 
n gi + ig 2 

8n /1 + if 2 
n gi + ig 2 



(13a) 
(13b) 



where 

di = {ax + a 2 + a 3 ) (AxC 2 - A 2 C X ) , (14a) 
d 2 = ai« 2 (C 2 - Cx) + axa^C2 1 (14b) 
fx = (a x +a 2 + a 3 ) (CxB 2 - C 2 B X ) , (14c) 
f 2 = —axa 3 C2 — a 2 a%Cx, (14d) 
.91 = -aia 2 a 3 + (ax + a 2 + a a ) (BxA 2 - AxB 2 ) , (14e) 
g 2 = axa 2 (Bx - B 2 ) + axa 3 (A 2 - B 2 ) + a 2 a 3 yli,(14f) 

with on = nuj/Xi (i — 1,2). The pressure oscillations can 
be given in terms of the instantaneous values of Sn, 8X ei 
and 8X S , 

dP 

8P=—Sn + n(C 1 -C2)8X e + nCxSX Sl (15) 
on 
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where the Ci are the same as in Eq. jSJ . In the derivation 
we took into account that rii = dPj d[ii and that the total 
pressure is given by the sum of the partial contributions 
of the quarks and electrons, P = ^ i Pi(ni). 

After taking into account the relation (fTS]) together 
with the solution for 5X e $ and 5X St o, see Eq. (fl3|) . the 
expression @ becomes 

(fidto) = ~Sn [(Cx - C 2 ) Im (<5X e , ) + Cilm (<SA S , )] . 

This shows that the imaginary parts of 5X e ^ and <Wf S) o 
are directly related to the dissipation of the energy den- 
sity. From the physical point of view, these imaginary 
parts are related to the phase shifts (i.e., lagging) in the 



oscillations of the electron and strangeness fractions with 
respect to the oscillations of the density of quark matter. 

Now, by comparing Eq. (|16p with the definition in 
Eq. ([1]) , we derive an explicit expression for the bulk vis- 
cosity, 

n 2 

C = -=— [(Ci - C 2 ) Im {SX ei0 ) + dim (SX S<0 )] . (17) 
ujouq 

Finally, after making use of the expressions for 5X e Q and 
SX Si0 , see Eq. (IIU)) , this result can be rewritten as follows: 

C = Ci + C2 + C 3 , (18) 

where 



[axa 2 a 3 C 2 + (ax + a 2 + a 3 ) (A Y C 2 - A 2 d) (A X C 2 - (Ax - B 1 )C 1 )] , 
[a ia2 a 3 Ci + ( ai + a 2 + a 3 ) [(A 2 - B 2 ) C x - (Ax - B X )C 2 \ \{A 2 - B 2 ) C x - A 2 C 2 ]] , 



n a 2 a 3 
^ ~ T~ — 2 

w 9i + g 2 

n a x a 3 

^ 2 — — r - ; — 2 

w 51 + 52 
n aiQ:2 

ft = q - ; — 2 

w 51 + 52 

This can be simplified using the relation A x ~ A 2 = B Xl which follows from the definitions in Eq 
obtain 



(19a) 
(19b) 



(1 



(2 



w 5? + 9% 
n a±a 3 

ugf + gl 1 
n a\a 2 

u 9i + 92 



a ia2 a 3 (Cx - C 2 y + (ax + a 2 + a 3 ) (BxC 2 - B 2 Cx) [(Ax - ^ 2 )C 2 - (B 2 -Bx + Ax- A 2 )Cx] (19c) 

. Therefore, we 

t3 r axa 2 a 3 C 2 + (ax + a 2 + a 3 ) (AxC 2 - A 2 Ci f , 

a x a 2 a 3 Cl + (ax + a 2 + a 3 ) \(A 2 - B 2 ) Cx - A 2 C 2 \ 
a x a 2 a 3 (Cx - C 2 ) 2 + (ax + a 2 + a 3 ) (BxC 2 - B 2 CxY 
I 



(20a) 
(20b) 
(20c) 



Let us emphasize that the Q, i = 1,2,3, are not di- 
rectly associated with the separate contributions of the 
three types of weak processes in Fig. [1] Because of an 
interference of the weak processes, each of these contri- 
butions depends on all three rates A^, i — 1,2,3. The 
separation occurs only in the high-frequency limit, i.e., 
ai = nio/Xi — > 00. Indeed, in this case 

d ~ ^irCf , (no sum over i), (21) 

where G 3 = C 2 — Cx- Note that the formal criterion for 
this separation reads u> S> Lo scp . An estimate for uj scp is 
derived from the parametric dependence of the coefficient 
functions As and B; on densities, 



y/Al (A 2 + A 3 ) 

2/3 



(22) 



In order to understand the general features of the re- 
sult in Eq. (|20|) . let us consider the interplay of the non- 
leptonic and semi-leptonic contributions to the bulk vis- 
cosity. In the limit of vanishing rates for the semi-leptonic 



processes, i.e., A 2 ,A3 — > 0, the bulk viscosity reduces to 
the following well-known result for strange quark matter 
[Sill: 



Cn 



AiC? 



lu 2 + (XxAx/nY 



(23) 



In this limit, there are only non-leptonic processes left, 
and they induce the dissipation of the oscillation energy. 

In order to better understand the interplay of different 
types of processes, it is instructive to consider also the 
limit of an infinitely large non-leptonic rate, Ai — > 00, 
keeping the semi-leptonic rates A 2 and A3 finite. In this 
case, the expression for the bulk viscosity is given by 



Clcp 



(A 2 + A 3 )(C 2 ~C 1 A 2 /^ 1 ) 2 
lu 2 + [(A 2 + A 3 ) (B 2 - BxA 2 /Ax) /n] 



(24) 



It may appear puzzling that the two seemingly equiva- 
lent limits, namely Ai 3> A 2 , A3, lead to such very differ- 
ent results. The problem can be resolved by noting that 
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the result in Eq. (|23[) is reliable only if the following ad- 
ditional constraint is satisfied: (Aa + A3) <C ti 4 ^ 3 lu 2 / 'Ai. 
[In deriving this constraint, we assumed that the coeffi- 
cient functions Ai and Bi scale with the density as n 1 / 3 , 
which will turn out to be a reasonable approximation, 
see Eq. (|27p below.] In contrast, the result in Eq. (|24]) 
is reliable only if Ai 3> n 4 / 3 uj 2 / (A2 + A3). It is conve- 
nient, therefore, to introduce the following characteristic 
frequency: 



LU 



yAj (A 2 + a 3 ) 

n a/3 



(25) 



which separates the two qualitatively different regimes. 
The result in Eq. (f23f is to be used when the frequency 
of the density pulsations is sufficiently high, uj 3> ujq, so 
that only the fast non-leptonic processes have a chance 
to dampen the kinetic energy efficiently. The result in 
Eq. (|24p. on the other hand, should be used at sufficiently 
low frequencies, w < wo, when only the much slower 



semi-leptonic processes provide a substantial damping of 
the oscillations. At intermediate frequencies, u> ~ loq, 
neither Eq. (|2"5|) nor Eq. (|24p provides a good approxima- 
tion to the bulk viscosity because of a strong interference 
of the non-leptonic and semi-leptonic processes. Below 
we study this interference in some more detail. 



III. BULK VISCOSITY IN NORMAL PHASE 

In order to calculate the bulk viscosity in the normal 
phase of three-flavor quark matter, we need to determine 
the corresponding thermodynamic coefficients A, , Bi and 
Ci [see Eq. ©] and calculate the difference of the rates 
of the three pairs of weak processes in Fig. [TJ 

Let us start from the derivation of the coefficients Ai, 
Bi and Cj. For this purpose, we use of the following 
zero-temperature expressions for the number densities of 
quarks and electrons: 



(/4 - mfj 



3/2 



2a s , 2 



m)) 



3m 2 



M 2 - m 2 



In- 



m f 



3vr 2 ' 



for / = u,d, s, (26a) 
(26b) 



Note that the expressions for quarks include the leading- 
order corrections due to strong interactions [H], [H, [H, 
[37j . By making use of these relations together with the 
definitions in Eq. I|8|). we derive the following approxi- 
mate expressions for the coefficient functions: 



Bi 
Cy 
A 2 
B, 

C-2 



2n 2 n 



3Md 



7r 2 n 








7r 2 n 




3~ 1 






mf. 


3^ s 


1 iL 

3/x u 



(27a) 
(27b) 
(27c) 
(27d) 
(27e) 
(27f) 



On the right-hand side of these expressions, we made 
use of the equilibrium relations satisfied by the chemical 
potentials, (j, s = /j,d = fx u + (i e , and neglected higher- 
order mass corrections in the expressions for Ai and Bi. 
In the expressions for Cj, on the other hand, we kept the 



leading-order mass corrections, which also include the 
contributions linear in the strong coupling constant a s , 



4a s m 2 



In 



2 M/ 



(28) 



3ttm/ V m f 

The reason for keeping the mass corrections is that the 
functions Cj vanish in the massless limit. In that limit, 
all £j's are zero as appropriate for scale invariant theories 

In order to further proceed with the calculation of the 
bulk viscosity in the normal phase of three-flavor quark 
matter, we also need to know the rate difference of the 
three pairs of weak processes in Fig. [TJ The rates of both 
the non-leptonic and semi-leptonic processes have been 
calculated in the literature, see Refs. [H, [H, H^, H(| and 
Refs. \^ , respectively. In the limit of three massless 
quarks, for example, the rates are 

&A tG 2 f cos 2 9 C sin 2 9 c ^ d T 2 , 



Ai 



A 2 



A3 



57T 3 

17 

40^ 
17 

15^ 



G% sin 2 9 c ^ s m 2 s T A , 
■G 2 F cos 2 9ca s HdHuHeT A . 



(29a) 
(29b) 
(29c) 



When the value of the strange quark mass is not small, 
the expressions for the rates are more complicated [H, Ho| 



and can only be calculated numerically. For the purposes 
of this paper it suffices to use the approximate expres- 
sions in Eq. (|29p. These provide a reasonable approxi- 
mation for the range of strange quark masses considered 
below. 

Now, by making use of the results for the rates (129)) 
as well as for the coefficient functions in Eq. (|2T|) . we 
can straightforwardly calculate the bulk viscosity. In the 
calculation, we use the following two representative sets 
of model parameters: 



set A set B 

n = 5po n = 10po 

m s = 300 MeV m s = 140 MeV 

a 3 =0.2 a s = 0.1 




In both cases the light quark masses are m u = 5 MeV 
and ma = 9 MeV. Regarding the choice of parameters 
for sets A and B, several comments are in order. Set 
B is supposed to be characteristic for the conditions in 
the inner core of a quark star. The strong coupling con- 
stant a s is small due to asymptotic freedom and the in- 
medium constituent strange quark mass m s assumes a 
value close to its current value on account of chiral sym- 
metry restoration at large baryon densities. Set A applies 
to intermediate densities, such as occur in the outer core 
of a quark star. Here, the strong coupling constant a s 
and the in-medium constituent strange quark mass m s 
assume larger values. 

It can be shown that the electron chemical potential is 
driven by both the strange quark mass and the leading- 
order corrections due to strong interactions, see Eq. (|26[) . 
Without a s corrections, the electron chemical potential 
would be 9.9 and 50.5 MeV for strange quark masses 
of 140 MeV (n = 10p ) and 300 MeV (n = 5p ), re- 
spectively. The leading-order corrections due to strong 
interaction tend to reduce the values of p e . A simple 
analysis shows that, in fact, the electron chemical poten- 
tial could even formally change the sign when the value of 
a s is sufficiently large (e.g., a s > 0.5). This would mean 
that strange quark matter requires the presence of the 
positrons rather than electrons to stay neutral. While 
not forbidden, such a possibility should be accepted with 
great caution. Indeed, the leading-order a s corrections 
may be unreliable in this regime. 

The general algorithm for calculating the bulk vis- 
cosity is as follows. First, by assuming a fixed value 
of the baryon density of neutral quark matter n (i.e., 
n = 5po ~ 0.75 fm or n = 10po ~ 1-5 fm -3 for the two 
cases considered), we determine the chemical potentials 
of the quarks and electrons in (3 equilibrium. For the two 
representative sets, the values of the chemical potentials 
as well as the coefficient functions Aj, Bi and are given 
in Table HI These are used in the calculation of the rates 
of the weak processes. Putting everything together, the 
result for the bulk viscosity follows from Eq. (|18p . 

The numerical results for the bulk viscosity as a func- 
tion of the period of density oscillations t are presented in 



i — i 




T[Sec] 



FIG. 2: (Color online) The bulk viscosity for the normal phase 
of three-flavor quark matter as a function of the period of 
density oscillations. Results for set B are shown in the upper 
panel and those for set A in the lower panel. For each tem- 
perature, the dots on the lines correspond to the values of the 
frequency defined in Eq. (1251) . 



Fig.[2]for the two cases: (i) n = 10po and m s = 140 MeV 
(upper panel), and (ii) n — 5po and m s = 300 MeV 
(lower panel). Different line types correspond to differ- 
ent values of the temperature: T — 0.1 MeV (solid lines), 
T = 0.2 MeV (dashed lines), T = 0.4 MeV (dotted lines), 
and T = 0.8 MeV (dashed-dotted lines). The thin lines 
show the high- and low-frequency approximations, de- 
fined in Eqs. (|2"3")l and (f2~4")l . for each value of the temper- 
ature. 

Note the double-step structure of the bulk viscosity 
as a function of r. The first step at small r is the usual 
one. It corresponds to the low-frequency saturation (with 

Cnolf^ ~ 1/Ai ) of the non-leptonic contribution to the 
bulk vicosity, see Eq. (I23|) . The second step at about 
t = 2tt/loq (marked by dots in the figures) is a quali- 
tatively new feature. As should be clear from our dis- 
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n/po 


m a [MeV] 


p e [MeV] 


p u [MeV] 


p d = p 3 [MeV] 


Ai [MeV] 


A 2 [MeV] 


Bi [MeV] 


B 2 [MeV] 


Ci [MeV] 


C 2 [MeV] 


5 


300 


39.139 


402.463 


441.602 


239.432 


127.937 


111.386 


-3.726 x 10 4 


-60.463 


-60.460 


10 


140 


7.396 


495.275 


502.671 


324.118 


164.288 


160.268 


-2.080 x 10 6 


-10.692 


-10.709 



TABLE I: Two sets of parameters used in the calculation of the bulk viscosity. 



cussion in the preceding section, its appearance is the 
consequence of the interplay between the non-leptonic 
and the semi-leptonic weak processes contributing to the 
bulk viscosity of strange quark matter. At sufficiently 
large r (i.e., low frequencies) the contribution of the semi- 
leptonic processes also saturates. 

The increase of the bulk viscosity at low frequencies 
due to the contributions of the slower semi-leptonic pro- 
cesses is not unexpected. (For the same reason, weak pro- 
cesses are much more important for compact stars than 
strong processes, operating on typical QCD time scales 
of order 1 fm/c). However, the main observation here 
is that the increase of the bulk viscosity due to the in- 
terplay between the non-leptonic and semi-leptonic weak 
processes could already be visible at frequencies relevant 
for the physics of compact stars. Moreover, by comparing 
the results in the two panels of Fig. [2], we find that for the 
conditions corresponding to lower densities, the range of 
frequencies where the semi-leptonic processes contribute 
widens significantly. 

The bulk viscosity as a function of the temperature is 
given in Fig.[3]for set B (upper panel) and for set A (lower 
panel) . The additional increase of the bulk viscosity due 
to the semi-leptonic processes is seen as a "bump" at in- 
termediate values of the temperature. The results in the 
two panels demonstrate once again that the relative role 
of the semi-leptonic processes increases with increasing 
period of density oscillations, provided the baryon den- 
sity is sufficiently small (set A). 

From Fig. O we find that the substantial modification 
of the bulk viscosity due to the Urea processes occurs at 
temperatures between about 0.1 MeV and 1 MeV. This 
temperature range is of relevance to young neutron stars 
and, therefore, should be studied in more detail. 

Perhaps the best way to appreciate the relative role 
of the Urea processes is to study the ratio between the 
complete expression for the bulk viscosity (fT5)) and the 
commonly used approximate form (|23p that takes only 
the non-leptonic interactions into account. In the more 
interesting case of parameter set A, this ratio is shown in 
Fig. 2] as a function of temperature for four different val- 
ues of the period of density oscillations: - = 1 Hz (solid 
line), i = 10 Hz (dashed line), ± = 100 Hz (dotted line), 
- = 1000 Hz (dashed-dotted line). The plot shows how 
the ratio C/Cnon changes when the period of oscillations 
t varies in the whole range from 1 ms to 1 s. 

The results presented in Fig. [3] are the main results of 
this paper. They show that neglecting the Urea processes 
can result in underestimating the value of the bulk viscos- 
ity by an order of magnitude, at intermediate densities 
and sufficiently large in- medium strange quark masses. 





FIG. 3: (Color online) The bulk viscosity for the normal phase 
of three-flavor quark matter as a function of the temperature 
for several fixed values of the frequency of density oscillations. 
Results for set B are shown in the upper panel and those for 
set A in the lower panel. 



We believe this finding might be of relevance for strange 
quark matter under conditions realized inside young neu- 
tron stars. 

Before concluding this section, it is appropriate to 
note that the role of the semi-leptonic processes is neg- 
ligible in the case of the smaller strange quark mass, 
m s = 140 MeV. This is seen from Fig. [3] (upper panel), 
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T [MeV] 

FIG. 4: (Color online) The ratio £/£ n on as a function of tem- 
perature for set A. The results for three fixed values of the 
density oscillation frequencies are shown: - = 1 Hz (solid 
line), i = 10 Hz (dashed line), ± = 100 Hz (dotted line), and 
~ = 1000 Hz (dashed-dotted line). 



and this can also be confirmed by studying the ratio 
C/Cnon- Ks value deviates from 1 by at most 28%. 
This is in qualitative a gree ment with earlier studies 
[H d, E3, IS GJ, HI M HI, neglecting the semi- 
leptonic processes. 



IV. DISCUSSION 

In this paper we have studied the subtle interplay be- 
tween the Urea (semi-leptonic) and the non-lcptonic weak 
processes in determining the bulk viscosity of neutral, (3- 
equilibrated strange quark matter. In general, the contri- 
butions of the two types of weak processes are not separa- 
ble. The exception is the high-frequency limit, uj 3> w scp , 
see Eq. (f22j) . in which case the contributions naturally 
separate. Because of a much higher rate of the non- 
lcptonic processes, they dominate in this high-frequency 
limit. 

With decreasing frequency the role of the Urea pro- 
cesses gradually increases. The maximal mixing with the 
non-leptonic processes occurs at about the frequency uiq, 
see Eq. (|25p for the definition. Depending on the specific 
choice of model parameters, typical values of ujq are in 
the range from about Is" 1 to 10 3 s~ x , see also Fig. [2] for 
the corresponding values of the period of oscillations. 

Our numerical results for the normal-conducting phase 
of strange quark matter demonstrate that the commonly 
used approximation, in which the Urea processes are 
completely neglected, could substantially underestimate 
the value of the bulk viscosity, see Figs. and [3] By 
sweeping over a wide range of model parameters, we find 
that the role of the Urea processes is most important in 
the range of temperatures between about 0.1 MeV and 
1 MeV. In the outer core of a quark star, i.e., when the 



density is not too large, the strange quark mass not too 
small, and/or the period of oscillations not too short, the 
inclusion of the Urea processes could lead to an increase 
of the viscosity by an order of magnitude, see Fig. |U 

The findings of this work could have important im- 
plications for the physics of young neutron stars with 
strange quark matter interiors and/or for pure strange 
stars that could potentially exist too. In connection to 
the r-modes instabilities driven by gravitational radia- 
tion, the increase of the bulk viscosity due to the Urea 
processes is likely to broaden the region of stellar sta- 
bility in the temperature-angular-frequency plane. The 
fact that the largest change occurs for the temperature 
range from about 0.1 MeV to 1 MeV might be also very 
important. Indeed, at the lower end of this range, the 
dominant role in suppressing the r-modes is expected to 
pass from the bulk viscosity to the shear viscosity. If 
the bulk viscosity is substantially higher than previously 
estimated, it may dominate the dissipative dynamics to 
much lower temperatures. 

In view of possible color superconductivity in strange 
quark matter, in the future one should also investigate 
the role of the Urea processes in color-superconducting 
phases with various types of spin-zero and spin-one 
Cooper pairing. Because of the gaps in the quasipar- 
ticle spectra, the interplay between the semi- and non- 
lcptonic processes in superconductors is expected to be- 
come more complicated due to the suppression of the 
rates. Certain general features of the corresponding dis- 
sipative dynamics could be predicted even without a de- 
tailed study. One could say, for example, that (i) the 
suppression of the rates should lead to the suppression of 
the bulk viscosity in the high-frequency limit, uj 3> w S cp, 
(ii) the border line itself, lu scp , should shift to a lower 
value, (iii) the rates of the Urea processes could possibly 
become even higher than the rates of the non-leptonic 
ones. The first two observations are simple consequences 
of our general results in Sec. |TTJ The situation described 
in item (iii) could be realized when all quark quasiparti- 
cles are gapped, for example, in one of the versions of the 
color-spin-locked phase. At low temperatures, the rates 
of the semi- and non-leptonic processes are suppressed 
by exponential factors exp(— <f>/T) and exp(— 2<p/T), re- 
spectively. A more detailed discussion of spin-one color- 
superconducting phases of strange quark matter is sub- 
ject of a subsequent paper (43j . 

Note added. While writing our paper, we learned that 
another study of the bulk viscosity of strange quark mat- 
ter is done by H. Dong, N. Su, and Q. Wang |4l |. 
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